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THE NEAR-IR CA II TRIPLET - a RELATION FOR BULGES OF SPIRAL GALAXIES 
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ABSTRACT 

We present measurements of the near-infrared Calcium II triplet (CaT, CaT*), Paschen (PaT) and Magnesium 
(Mg I) indices for a well-studied sample of 19 bulges of early to intermediate spiral galaxies. We find that both the 
CaT* and CaT indices decrease with central velocity dispersion a with small scatter. This dependence is similar to 
that recently found by Cenarro (2002) for elliptical galaxies, implying an uniform CaT* - a relation that applies 
to galaxies from ellipticals to intermediate-type spirals. The decrease of CaT and CaT* with a contrasts with 
the well-known increase of another a-element index, Mg2, with a. We discuss the role of Ca underabundance 
([Ca/Fe]<0) and IMF variations in the onset of the observed relations. 

Subject headings: galaxies: bulges — galaxies: fundamental parameters — galaxies: abundances 



1. INTRODUCTION 

Scaling relations provide powerful insights on the forma- 
tion and evolution of early-type galaxies. Besides linking in- 
ternal structure and dynamics through the Fundamental Plane 
(Djorgovski & Davis 1987), they allow to link these proper- 
ties with those of stellar populations. The most studied scaling 
relation between population and dynamics is the Mg2 - a re- 
lation (Terlevich et al. 1981), which shows that the Mg2 line 
index increases with galaxy mass. This scaling relation may 
be used to study whether bulge populations are younger than 
those of elliptical galaxies, with the basic assumption that el- 
liptical galaxies in clusters are old. However, stellar popula- 
tion studies in the wavelength region of the Mg2 index (Trager 
et al. 2000; Kuntschner 2000) indicate that Mg 2 is affected by 
age and metallicity; that Mg is overabundant w.r.t. Fe in bright 
galaxies; and, therefore, that the Mg2-er relation for ellipticals 
is not necessarily a relation for old galaxies (Trager et al. 2000). 

The near-IR Ca II triplet at 8500 A is a bright atomic fea- 
ture that can easily be studied with modern detectors. It has 
been used in the past to analyse stellar populations (e.g. Carter 
et al. 1986; Terlevich et al. 1990; Olszewski et al. 1991; Bica 
et al. 1998). Recently, a series of papers have been published 
describing a new stellar library that is used to make a thorough 
calibration of the Ca II triplet and other features in the same 
spectral region in terms of stellar parameters (Cenarro et al. 
2001a,b, 2002). A followup paper (Vazdekis et al. 2003, here- 
after VAZ03) presents the behaviour of these features for single 
stellar populations as predicted by recent evolutionary synthesis 
models. Cenarro (2002, hereafter CEN02), have used this cali- 
bration to study the near-IR Ca II triplet in a sample of elliptical 
galaxies. They find that CaT*, the Ca triplet index corrected 
for the presence of Paschen lines, decreases with a, and that, in 
the most luminous galaxies, CaT* is much lower than predicted 
by models. The decrease of CaT* with a may be a result of 
a dwarf-enhanced IMF or a low Ca/Fe abundance ratio for gi- 
ant ellipticals. In a similar study, Saglia et al. (2002, hereafter 
SAG02) find a similar underabundance of CaT w.r.t. models 
in luminous ellipticals, but do not see any dependence of CaT* 
with a. 

1 also CRAL, Observatoire de Lyon, F-69561 St-Genis Laval cedex, France 



Here we analyse a sample of galactic bulges to compare the 
trends of Ca with a in bulges to those found in ellipticals. Dif- 
ferent abundance ratios might be expected if bulges and ellipti- 
cals have different formation histories. The study offers an op- 
portunity to see whether CaT* - a relation may be established 
in general for the spheroidal components of galaxies. 

2. OBSERVATIONS & DATA REDUCTION 

Nineteen early- to intermediate-type galaxies (SO to Sbc), 
from the sample defined by Balcells & Peletier (1994), were ob- 
served using the ISIS double spectrograph at the 4.2m William 
Herschel Telescope of the Observatorio del Roque de los 
Muchachos at La Palma in 1997 July 11-13. Minor-axis spec- 
tra were obtained in the 8360 - 9170 A range, at 1.74 A spec- 
tral resolution (FWHM of the arc lines). The slit width was 
1.2" to match the seeing at the time of the observations. A de- 
tailed summary of the observations and data reduction can be 
found in Falcon-Barroso et al. (2002). We measured the re- 
cently defined line-strength indices CaT*, CaT, PaT (Cenarro 
et al. 2001a) and Mgl (CEN02), on the flux calibrated spectra. 
The main advantages of these new definitions over the previous 
ones (e.g. Armandroff & Zinn 1988; Diaz et al. 1989; Rutledge 
et al. 1997) are that they have been designed to ease problems 
in the definition of the continuum bands, the contamination of 
the Paschen lines, and the sensitivity to the S/N ratio and ve- 
locity dispersion (see Cenarro et al. 2001a for more details). In 
order to explore the dependence of the relation of line-strength 
and velocity dispersion on aperture, we measured the 4 line- 
strength indices in three different apertures (r e ff , r e ff/2, 1 .2 x 4 
arcsec 2 ). The signal-to-noise ratios for the different apertures 
range from 67 to 256, 55 to 232 and 75 to 244 per angstrom, 
respectively. if-band effective radii were taken from Table 1 
of Falcon-Barroso et al. (2002) and converted to the minor axis 
using high-spatial resolution iJ-band ellipticity profiles from 
Balcells et al. (2003). For a number of our bulges it was not 
possible to accurately measure the index at aperture sizes < 
r c ff/2 due to its size and seeing conditions. The difference at 
smaller apertures, however, should be minor since the near-IR 
Ca II triplet depends very little with the size of the aperture. For 
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brevity, only the r e fj, and 1.2 x 4 apertures are shown in Fig. 1. 
The measured values for the r c ff/2 aperture have been included 
in the electronic version of Table 2. Kinematic parameters were 
derived using FOURFIT developed by van der Marel & Franx 
(1993). A gaussian fit to the LOSVD was performed to ob- 
tain mean radial velocities and velocity dispersions. We have 
made use of a single-age, single-metallicity synthetic model 
([Fe/H]=0.0, t=6.31 Gyr) from VAZ03 as a kinematic template. 
This procedure has proven to be very effective in reducing tem- 
plate mismatch and produces as good results as the best stel- 
lar templates in this wavelength domain (Falcon-Barroso et al. 
2003). Tabulated values for the CaT* index and velocity dis- 
persions along the minor axis of the galaxies can be found in 
the same paper. We broadened all the spectra to 300 km s _1 , 
slightly higher than the largest velocity dispersion in our sample 
(291 km s" 1 for the central aperture of NGC 5838). This cor- 
rection is small for the CaT* index, which has shown to be fairly 
constant with velocity dispersion up to 300 km s _1 (VAZ03). 
Errors in the indices include contributions from radial velocity 
uncertainties and Poisson errors (see Falcon-Barroso et al. 2003 
for more details). The uncertainty from S/N dominates the er- 
rors in the indices. We have corrected our indices for differ- 
ences between the spectro-photometric system of the data and 
the model predictions using a set of 10 stars in common with 
the stellar library on which the models are based. The following 
offsets are found (Models-Us): ACaT* =0.178±0.032, ACaT 
=0.205 ±0.027, AMg I =-0.001 ±0.022, APaT=0.029±0.007. 



Table 1 
Fitting Parameters 



Index 


a 


Aa 


b 


Ab 


r 2 


r s 


a(r s ) 


rms 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 










1.2x4 










CaT* 


8.72 


0.38 


-0.80 


0.17 


0.51 


-0.71 


7.4E-04 


0.09 


CaT 


8.89 


0.40 


-0.70 


0.19 


0.41 


-0.58 


8.6E-03 


0.07 


Mgl 


0.00 


0.08 


0.11 


0.04 


0.13 


0.20 


4.2E-01 


0.03 










r c ff 










CaT* 


9.44 


0.39 


-1.14 


0.18 


0.52 


-0.76 


1.7E-04 


0.11 


CaT 


9.67 


0.42 


-1.06 


0.20 


0.63 


-0.80 


4.2E-05 


0.07 


Mgl 


0.04 


0.09 


0.10 


0.04 


0.13 


0.22 


3.6E-01 


0.02 



NOTE: Parameters of the error-weighted least-squares linear fits to the index 
- log(cr) relations, defined as I = a + b- log(cr) at 300 km s _1 resolution. 
Columns 2, 3, 4, 5 show the fitted parameters and their errors. Columns 6, 
7, 8 give the linear correlation coefficients (r), the Spearman rank-order cor- 
relation coefficients (r s ) and significance levels (a(r s )). Column 9 lists the 
rms around the best fit. 

3. THE INDEX - LOG(ct) RELATIONS 

The relations of the line-strength indices CaT*, CaT, PaT, 
Mg I and velocity dispersion are shown in Figure 1 . We present 
each relation for the two aperture sizes described in § 2 (labeled 
on the top boxes). Early-type spirals (T < 0) are depicted with 
filled symbols, and later types (T >0) with open symbols. Index 
values for a range of SSP (single-age, single-metallicity stellar 
populations) models at the resolution of the data are shown in 
the fourth-column panels. Linear least-squares fits are plotted 
for some of the relations, with best-fit coefficients given in Ta- 
ble 1 . The data values are presented in Table 2. 

A strong anti-correlation exists between the CaT and CaT* 
indices and log(er) in all apertures. This is reflected in the low 
value of the significance levels of the Spearman rank-order tests 
(Table 1, column 8), which give the probability of a null corre- 
lation, and the rms values about the best fit for the three indices. 
The anti-correlation of the CaT* index with a shown here for 
bulges is similar to that recently found by CEN02 for ellipti- 
cal galaxies (which was, if at all, marginally detected also for 



ellipticals by SAG02). 
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FIG. 1 . — The CaT* , CaT, Mgl and PaT index relations with log(cr) for dif- 
ferent apertures. The data has been convolved to a spectral resolution of 300 
km s . Solid circles are early-type bulges (T < 0), while later-type bulges 
(T>0) are represented by open circles. Solid lines represent an error-weighted 
least-squares fit to the all the points. See Table 1 for more details. The last 
column shows SSP model predictions for different ages and metallicities, for 
a Salpeter IMF, for each index. Dashed, solid, dotted and dash-dotted lines in 
the last column represent SSP models of [M/H]=0.2,0.0,-0.4,-0.7 metallicities 
respectively. The thick solid line, in the last column, shows the effect of a 
higher IMF (i.e. 2.3) for a [M/H]=0.0 model. 

The Call triplet is the first metal line index that is found 
to anticorrelate with a. The Mgl index, in contrast, shows 
a positive trend with cr as expected from the behaviour of the 
Mgfe and Mg 2 indices in the optical range (e.g. Terlevich et al. 
1981). We have investigated whether different apertures would 
lead to significant differences in the index-cr relations. No ma- 
jor difference is found between the gradients of any index at 
each aperture: the slopes and zero-points agree within the er- 
rors (Table 1). It is interesting to see that none of our bulges 
give a measurement of the CaT or CaT* indices that lies above 
the predictions of the models (VAZ03). This result agrees with 
the model predictions that the CaT and CaT* indices tend to 
saturate for metallicities above [M/H]— 0.4 (VAZ03), which 
are also supported by the data of CEN02 and SAG02. Given 
that our sample contains bulges from SO to Sbc, we can study 
whether the index-cr relations depend on morphological type 
(see Fig. 1). Due to our limited number of galaxies, we split 
the sample in just two groups: one of them containing galax- 
ies of types SO, SO/a (T<0) and a second group with the Sa to 
Sbc galaxies (T>0). No differences are found for any of the 
relations. 

We have compared our index-cr relations with those found 
by CEN02 measured for an aperture of 2x4 arcsec 2 . An accu- 
rate comparison with SAG02 is not possible since their aper- 
ture is four times smaller than the smallest we can get given 
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the size of some of our bulges. We perform the comparison 
with CEN02 at a similar aperture size and at the same spec- 
tral resolution (400 km s _1 ). In figure 2 we show the com- 
parison of the CaT* index. The results are in good agree- 
ment at the 1-cr level, with the slopes agreeing within the errors 
(b[Us]=-0.69±0.18,b[CEN02]=-0.56±0.15). The exclusion of 
one of our low-velocity dispersion galaxies (NGC 7457), pro- 
vides an even better agreement in the slope values (b[Us]=- 
0.52±0.19). However the rms of our data around the best fit 
is much smaller (rms[CaT*(Us)]=0.14) than the one measured 
by CEN02 (rms[CaT*(CEN02)=0.28). On the other hand, our 
sample shows values for the PaT index slightly higher than 
those of CEN02 for ellipticals (-0.03 A at 400 km s" 1 reso- 
lution). Previous studies of the same sample by some of us, 
studying deviations of bulges from the FP of ellipticals and SOs 
(Falcon-Barroso et al. 2002) or color-color diagrams (Peletier 
et al. 1999a), conclude that bulge populations are old, but could 
be younger than ellipticals by up to 2.5 Gyr. The higher PaT 
values, and correspondingly younger ages for bulges could also 
be due to young stars in the centre of our galaxies. However, 
model predictions indicate that young populations (<1 Gyr), 
which increase the value of the PaT index, also lower the value 
of CaT*, which is opposite to what we find here. Given the 
small significance of the offset it is more likely that flux cal- 
ibration differences contribute to the observed shift. For the 
most extreme case in the sample (NGC 5965,PaT=0.91 A) the 
high PaT could be due to star formation induced by the pres- 
ence of a strong bar in this galaxy (Kuijken & Merrifield 1995). 
On the other hand, we find no correlation of the PaT index 
with (j, a result supported by CEN02 for Es, and contrary to 
the strong trend found by SAG02. An explanation for that dif- 
ference could be in the velocity dispersion correction they have 
applied to each galaxy in their sample. A large correction for 
more massive galaxies could explain the positive correlation of 
the PaT index and the flatness of the CaT* index as a function 
of velocity dispersion. 



Table 2 
The Data 



Galaxy 


log(o-) 


CaT* 


CaT 


Mgl 


PaT 


Aperture 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGC 5326 


2.224 


6.897 


7.238 


0.252 


0.367 


r e fF 




(0.008) 


(0.116) 


(0.124) 


(0.026) 


(0.105) 




NGC 5389 


2.091 


6.748 


7.363 


0.260 


0.661 


r c ff 




(0.010) 


(0.132) 


(0.141) 


(0.029) 


(0.119) 




NGC 5422 


2.221 


6.559 


6.951 


0.286 


0.422 


r c ff 




(0.010) 


(0.117) 


(0.125) 


(0.026) 


(0.106) 




NGC 5443 


2.002 


7.016 


7.227 


0.224 


0.228 


r c ff 




(0.015) 


(0.208) 


(0.223) 


(0.046) 


(0.189) 





NOTE: The log(cr), CaT*, CaT, Mgl and PaT values for for our sample of 
bulges in the different apertures, shown in Fig. 1. Uncertainties in the pa- 
rameters are shown in brackets. The following offsets must be added to the 
tabulated values in order to obtain the values shown in Fig. 1 : ACaT* =0. 178, 
ACaT=0.2()5, AMg I=— 0.001, APaT=0.029. The complete version of this 
table is in the electronic edition of the Journal. The printed edition contains 
only a sample. 

4. DISCUSSION 

The origin and abundance evolution of the atomic element Ca 
has been a matter of debate. Due to its nature as an a element, 
an evolution history similar to that of Mg is predicted (Woosley 
& Weaver 1995). However several observations reveal a more 
complex picture where Ca and Mg originate and evolve differ- 
ently. The problem has resuscitated again recently in the light 
of CEN02, SAG02. For early-type galaxies, the Ca abundance 
does not follow Mg and there is evidence that it is not enhanced 



to Fe, from, e.g. measurements of the Ca4227 line (Vazdekis 
et al. 1997, 2001) or the near-IR Ca II triplet (CEN02; VAZ03). 
On the other hand, the Mg - log(cr) relation has been interpreted 
as a mass-metallicity relation for early-type galaxies (Terlevich 
et al. 1981). In the light of the new results obtained by Cenarro 
et al. (2003) and model predictions (VAZ03), such statement 
does not necessarily hold for the near-IR Ca II triplet - log(er) 
relation (see Worthey 1998 and Henry & Worthey 1999 for ex- 
tensive reviews). 

The large discrepancy between the observed CaT* and val- 
ues for a Salpeter-IMF model for galaxies with large velocity 
dispersions can be interpreted as underabundance of Ca w.r.t. 
the other metals in more massive galaxies. Other possibilities 
are: 1) An IMF biased towards low-mass stars, 2) Composite 
Stellar Populations including a low metallicity component or a 
young and old population of the same metallicity, and 3) Ca de- 
pletion. We discuss the advantages/disadvantages of each one 
of these possible explanations in turn. 

The problem of the Ca underabundance has been previously 
addressed by several authors (O'Connell 1976; McWilliam & 
Rich 1994; Vazdekis et al. 1997; Peletier et al. 1999b; Vazdekis 
et al. 2001, VAZ03). It is hard to find an evolutionary scenario 
that can produce non-solar Ca/Mg values. Assuming that both 
Ca and Mg are produced by two "flavours" of SNe in which 
more massive stars (20-40 M Q ) produce Mg before («2.6 Myr) 
Ca, generated in less massive stars (12-30 M©), one could ex- 
plain a small discrepancy in abundances between Mg and Ca 
(Molla & Garcia- Vargas 2000). However, large abundance ra- 
tio differences are hard to obtain in this way, and impossible 
to get for real galaxies, for which the formation timescale is 
probably much larger than a few Myr. If we would use a- 
enhanced isochrones (Salasnich et al. 2000; VandenBerg et al. 
2000,VAZ03), required to get the correct [Mg/Fe] ratios in gi- 
ant ellipticals, assuming a solar-like element partition for the 
[Mg/Ca] and [Fe/Ca], we would predict even larger Ca indices 
(VAZ03), increasing the discrepancy. 

Next we discuss the apparent Ca underabundance at large ve- 
locity dispersions in terms of a dwarf-enhanced IMF. The uni- 
versality of the IMF has been a matter of controversy in the 
last decades (see Gilmore 2001 and Eisenhauer 2001 for cur- 
rent reviews of observational keys in favour and against this 
statement). IMF time variations among ellipticals have been 
studied by Cenarro et al. (2003) in detail. An advantage of 
this approach is its ability to explain both, the large Mg and 
the [Mg/Fe] values observed in giant ellipticals (Vazdekis et al. 
1996; Cenarro et al. 2003). The effect of a dwarf-dominated 
IMF is a weakening of the CaT* index (VAZ03; Cenarro et al. 
2003). This effect is stronger at high metallicities. As also 
mentioned in SAG02, there are four observational drawbacks 
against this hypothesis: first: indications of an IMF flatter than 
Salpeter in the bulge of our own Galaxy (Zoccali et al. 2000), 
second: the fact that the predicted FeH 9916 A feature is sig- 
nificantly stronger than what it is observed (Carter et al. 1986), 
and third: the discrepancy between visual M/L ratios (M/L> 10, 
Vazdekis et al. 1996; Maraston 1998) and dynamical estimates 
M/L«6 (Gerhard et al. 2001). An additional drawback is that 
visual-infrared (V — K) colors are just at the verge of what it is 
observed in Es, as pointed out by Cenarro et al. (2003). 

The use of composite stellar populations is analysed by 
SAG02. Their analysis reveals that a combination of popula- 
tions including a low-metallicity component was not able to 
explain the low value of CaT* and CaT together with H/3 for 
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the most massive galaxies. Our independent analysis confirms 
their result. 

Ca is known to be highly depleted onto dust in the interstel- 
lar medium (Crinklaw et al. 1994). Could it be that it is also 
depleted in stars? The evidence is not strong. Such depletion 
in stars would require that, during star formation, interstellar 
dust is recycled less than interstellar gas. Moreover, depletion 
would have to be stronger in more dusty galaxies, so one would 
expect a large scatter in the CaT* - a relation due to depletion. 
At least, the scatter for bulges of spirals, which are generally 
more dusty than ellipticals, should be larger. Since this is not 
seen, we argue that depletion can be excluded as the origin of 
the low CaT* values. 

In summary, many mechanisms can be used to try to explain 
the apparent underabundance of Ca for massive galaxies. How- 
ever none of them alone seems to be able to explain all the 
observations. In the light of the discrepancies found between 
all the possibilities, a combination of several effects could be in 
place. More work is still needed to provide an answer to this 
problem. 

We thank Javier Cenarro and Javier Gorgas for providing us 
with software and data in electronic form, and for many use- 
ful discussions that have helped to improve this letter. JFB ac- 
knowledges the finantial support from a PPARC studentship. 
The William Herschel Telescope is operated on the island of La 
Palma by the Isaac Newton Group in the Spanish Observatorio 
del Roque de los Muchachos of the Instituto de Astroffsica de 
Canarias. 
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FIG. 2. — Comparison of the CaT*- log(cr) relation between our sample 
of bulges and early-type galaxies from CEN02. Solid circles and line repre- 
sent our sample of galaxies and an error-weighted least-squares linear fit to it. 
Open circles and dotted line represent the CEN02 sample of ellipticals and its 
fit. The dashed line is an error-weighted least-squares linear fit to our sample 
excluding the galaxy marked with a circle (~1.8 in log(cr)). In order to com- 
pare our results we have degraded the resolution of our spectra to 400 km s ~ 1 . 
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